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abstract ; 


Following the big -bang baryon symmetric cosmology of Omnes 
where an initial phase separation of matter and antimatter leads 
to regions of pure matter and pure antimatter containing masses of 
the size of galaxy clusters by a redshift which we calculate to be 
of the order of 500-600, we show that at these redshifts , annihilation 
pressure at the boundary between the regions of matter and antimatter 
drives large scale supersonic turbulence which can trigger galaxy form- 
ation, This picture is consistent with the y-ray background observations 
discussed previously by Stecker, Morgan, and Bredekamp. Gravitational 
binding of galaxies then occurs at a redshift of ^ 70 at which time 
vortical turbulent velocities of ~ 3 X 10^ cm/s lead to angular momenta 
for galaxies comparable with measured values. 
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INTRODUCTION: 

In this paper, we will attempt to present a general schame which we feel to 
be a plausible model for galaxy formation deserving of future study. This scheme 
draws upon and attempts to synthesize three basic concepts relating to big -bang 
cosmology and galaxy formation. These concepts are:(l) that galaxy formation may 
have been triggered by turbulence in the cosmic medium at an early stage in the 
evolution of the universe, (2) that the universe in its early stages was a& emulsion 
of regions containing essentially pure matter and pure antimatter, and (3) that 
the observed background flux of cosmic gamma -radiation can be explained as due to 
matter-antimatter annihilation taking place on the boundaries of these regions 
with the observations providing a measure of the annihilation rate having occurred 
in the distant past. The first concept was examined by von Weizsifcker ( 1951) and 
gained futher importance when Gamow pointed out that large density fluctuations 
lik& those produced by turbulence were needed to trigger galaxy formation in an 
expanding universe. Gamow (1954) then further developed this argument and also 
pointed out that fossilized turbulence could account for the nonuniform spatial 
distribution of galaxies. Oort further pointed out that a large scale turbulent 
eddy might bring matter together in a manner that would result in associations 
of galaxies in systems with a positive total energy (Oort 1969) . Further contribu- 
tions to theories of cosmological turbulence were made by Nariai (1965a, b), 

Ozernoi and Chernin (1968), Sato, et al.(1970), Peebles (1970), Ozernoi and 
Chibisov (1970) and Harrison (1970a, b , 1971) . In all of these previous discussions, 
the initial turbulence had to be introduced in an ^ hoc manner; we will show here 
a natural way of feeding turbulent energy into the cosmic medium as a result of 
antimatter annihilations following concept (2) . 

The second concept was suggested by Harrison (1967) to reconcile big -bang 
cosmology with the principle of baryon symmetry, but was introduced by him as 
an ^ hoc initial condition. However, it was later reintroduced by Omnes in a 
much more satisfactory manner (Omnes 1969). Omnes suggested on the basis of strong- 
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interaction theory, that if we examine the big -bang model in its very early high- 
density state (corresponding to a critical temperature of ^50 MeV) a two -component 
phase separation could occur between matter and antimatter on a macroscopic scale. 

In subsequent papers (Omn&s . 1970, 1971a, b) . he traced the subsequent evolution 
and growth of the suggested matter -antimatter emulsion and showed that regions 
of pure matter and pure antimatter in the picture could grow to the size of galaxies 
or clusters of galaxies. 

The third concept was examined by Stecker, et al. (1971 ) with the Harr is on - 

Omnes picture in mind . They showed that such cosmological models could lead to the 
production of a background spectrijm of cosmic gamma -radiation having the same features 
which seem to be present in the observed spectrum. 

We will attempt to show here that the Omnes model which leads to annihilation 
rates of the order of those compatible with the gamma -ray observations also leads 
to annihilation pressures (at an earlier stage) which are sufficient to produce 
large scale turbulence with a maximiiim eddy size of the order of g^axy clusters. 

We then apply the general formulation of Ozernoi and his collaborators (but re- 
placing his ad hoc "photon -eddy” hypothesis by annihilation-generated turbulence) 
to show how annihilation -generated turbulence can trigger galaxy formation and 
lead to angular momenta for galaxies comparible with measured values. 
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II: THE OMNES MODEL: 

Omnes , following Gamov (1948a, b), considers a big-bang model which is initially 

at a very high temperature and density ( T oc r ^ where T is the cosmic temperature 

and R is the geometric scale factor of the universe ) . Bosons and fermion pairs 

exist in statistical equilibrium with photons of the blackbody radiation field at 

2 

a temperature T Me /k where M is the rest -mass of the boson or fermion considered. 
Clearly, at a high enough temperature, the number of bar yon -antibar yon pairs in the 
universe is of the same order as the number of photons in tfermal equilibrium. However, 
at present (T?:s2.7 K) , the number of baryons is a small fraction of the number of. photons 

7] = N„/N ^ 10"^° - 10“® 

D Y 

but this fraction is still much larger than that which would remain in thermal equil- 

-8 

ibrium in a blackbody at 2.7K (Chiu 1966). In fact, the ratio T) = 10 is reached in 

a blackbody at T 30 MeV. The two possible ways to account for this discrepancy are 

-8 

1) That there was initially a small excess of baryons of the order of 10 of the 
antibaryon number. 2)That the universe is baryon-S 3 ncnmetric and that a separation of 
matter from antimatter occurred at T >30 MeV. A consistent theory leading to the 
fulfillment of the second hypothesis based on particle physics and statistical argu- 
ments has been given by Omnes (1969) and further refined by Aldrovandi and Caser 
(1971). In this model, various mesons are assumed to be bound states of the N-if 
system with the appropriate quantum numbers. In a ** lattice gas” model then, a nucleon 
and an antinuefeon when placed within a cell of size r -- 1 fermi, look like a meson. 
However, the total number of meson is determined by statistical equilibrium at temp- 
erature T so that in a statistical average nucleons and antinucleons are excluded 
from occupying the same cell together. No such statistical exclusion principle applies 
to like baryons. Such a gas is similar to the model presented by Widom and Rowlinson 
{1910y for the study of liquid-vapor phase transitions. Their model is thermodynamically 
equivalent to a two-component system in which the pair -potential between molecules of 
like species is zero while that between unlike species implies a mutually excluded 
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3 

volume ^ . In such a model, Ruelle (1971) has rigorously shown the existence of 

a phase separation between the like and unlike species. At a stage when the baryonic 

-3 

density of the universe n ^ r or more (corresponding in Omnes” terminology to 

B 

•k 

T ^ ) this phase separation can then be expected to take place . Daschen, et al. 

(1969) have generalized the formalism originally given by Beth and Uhlenbeck (1937) 
to express the second virial expansion coefficient of the free energy of a system 
of strongly interacting particles in terms of observed scattering phase shifts. 

Omnes (1970) applied this treatment in his original work to show that above a crit- 
ical density, the free energy of a gas of nucleons and antimnucleons is a minimum 
for a phase of separately existing .components containing excesses of nucleons and 
antinucleons respectively. The work of Ruelle proves that such a phase separation 
will indeed occur. 

Following the epoch in the expansion of the universe when the density drops 
below the critical value (T<T^) , the two phases become unstable and matter and 
antimatter tend to mix and annihilate. Nevertheless, this process is slowed down 
by the strong pressure produced by annihilation at the boundary between regions 
of opposite baryonic number (the Leidenfrost effect pointed out by Alfven and 
Klein (1962)). During this annihilation period, the characteristic size of the 
regions of excess baryons and excess antibaryons in the emulsion is roughly equal 
to the diffusion length of the baryons. Omnes defines the end of the "annihilation 
period as corresponding to the time when the radiation temperature drops to ~ 30 
keV . At this point enough of the electron -^positron pairs in the blackbody rediation 
have annihilated so that the mean-free -path of high energy y-rays (produced by the 
This statistical argument, which is different from the one given by Omnes in the 
references, was pointed out by him to one of us in a private discussion. 
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decay of rr^-mesons arising from baryon-antibaryon annihilation) increases to the 
size of an average emulsion cell. At this point then, the annihilation pressure 
gradient extends over the whole cell and momentum can be transmitted on macroscopic 
scales of the order of the cell size to the fluid as a whole. Velocities can then 
be induced as large as 


V (nM c^/e)v. (1) 

where n is the baryon density, e is the blackbody radiation density and v^ is 
the thermal velocity of the baryons in temperature equilibrium with the black- 
body radiation ( Omnes 1971c). From 

this point on, in what Omnes defines as the coalescence period (Omnes 1971b). 
the characteristic size of the emulsion regions, d, grows in a way such that 

d ^ X (2) 

Y 

where \ is the mean free path of the p ion -decay v-rays . Equation (2) holds 
Y 

until the point when the radiation energy density (e) drops to the value equal 

2 - 8 /3 -2 

to that of the matter energy density (nM^c ) (Note: e oc t and n oc t at 

this time, where t is the time after the initial big-bang) . This corresponds 

4 

to a time 1.5 X 10 yr . after the big -bang when the corresponding blackbody 

4 

temperature was ~ 3.2 X 10 K (aee Appendix 1 for cosmological relations). 

The important relation giving the size of the emulsion regions of matter and 
antimatter as a function of blackbody temperature is shown in Figure 1, For 
further details of the Omnes cosmological model, we refer the reader to the 
papers by Omnes listed in the references. 
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III; THE PIASMA PERIOD 


(~5 X 10^ 


< z ^ 


10 


) 


We begin our original discussion of galaxy formation where Omnes'* 

4 

coelescence model stops, i.e., -- 1.5 X 10 yr after the ”big-bang” at a 

4 

redshift z 1.2 X 10 when the matter density p is equal to the radiation 

2 4 

density e. The following period ( ~ 5 X 10 ^ z < ~ 10 ) when p > e, but 

when the cosmic matter is still in the plasma state, we designate the plasma 
period (see appendices -1 and 2). 

The study of annihilation on a boundary region in appendix 2, shows that 

in a macroscopic collision of matter and antimatter regions with velocity 

the matter contained within a layer of thickness X from the boundary will 

X 

move away from the boundary with a velocity 

(2Fcv^P* (3) 


where F is the d 3 mamic efficiency factor given by equations (A2-30) and 
(A2-31). It follows from equations (A2-30) , (A2 -3 1) , and (3) that if we con- 
sider an initial estimate of v^^ to be roughly given by equation (1) for the 

end of the radiation period (z=z ) so that v.^^ v^ 2 X 10 cm/s, v_^ > v-. 

eq fl t ’ f2 fl 

and the collision is superelastic. Referring to appendix 2, if the total 
mass of plasma, , in a fluid cell involved in a collision is larger than 
M- , the mass of the cell within a distance X of the boundary, v_^ will in 
general be less than the value given by equation (3). The average fluid vel- 
ocity will increase in boundary collisions only if 


d/\. 


^X 




(4) 


estimated by a simple plane geometry model. 


When the plasma moves with a random fluid velocity v^ ^ ^exp 
where is the cosmological expansion velocity, the average coelescence 

growth rate for emulsion regions (dd/dt) v^ because mat ter -antimatter 
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collisions are elastic or superelastic and matter -matter collisions are in- 
elastic. Thus, at a time t(z), the characteristic emulsion' -size is of the 
order 

d(z) ~ v^(z)t(z) (5) 

where t(z) is given by (Al-7). 

At the beginning of the plasma period, the annihilation rate must be 
at least equal to the value given by the diffusion process and this value 
for is given by equation (A2-7). The average annihilation rate per unit 
volume is then 

Y - Y /d (6) 

V s 

Combining (6) and (A2-7) we find 

Y =-0.37N^v^/d (7) 

V 0 t 

At z - z , V as given by equation (1) is smaller than dX /dt 
eq max o ^ ^ 

(the coelescence growth rate during the radiation period as given by Omn^s) 

by almost a factor of 30, but the contributions of the Omnes coelescence 
process and the collision process (appendix 2) in producing macroscopic 
motions must at least be of the same order so that 

^ 2 (8) 

The value of in equation (7) is ^ n/20 (Omnes 197lt>)so that the 
lower limit on Y^ given by pure diffusion is 

Y ,=“3.3 X 10"^^(l+z)^'^ 

(9) 

for z z 

eq 

The maximum dynamic fluid velocity produced by annihilation is given by 

= [(0.27 -2mpC^Y^^^Ft)/(| iiipn)]S 
1 

==- c( f ,Ft/n ) 2 
v,d 


( 10 ) 
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so that for z 2 X 10 cm/s. This value is much higher than and 

we conclude that the annihilation rate Y will be determined by fluid motions 
rather than diffusion (at these redshifts) and ^ later stages of 

the plasma period, most of the annihilation results from fluid collisions with 
a velocity of the order of the average fluid velocity and we thus expect the 
minimum fluid velocity to be given by the relation (10). Using condition (4) 
together with equation (3) at z ~ 600 when F ~ 1, we find the upper limit 
for d at the end of the plasma period to be 

“„.x ~ 


The corresponding fluid velocity as given by equation (5) is 


13 - 3/2 

V- ^2,5 X 10 (1+z) cm/s 

f ,max 


( 12 ) 


A minimum value for d at the end of the plasma period can be obtained 
from equations (5) and ( 10) ; 

d . =-2.3X (13) 

min 


When d reaches the limit d given by equation (11), v cannot be in- 

max r 

creased further so that by combining equations (10) and (12) we find for an 
upper limit on the annihilation rate at the end of the plasma period when F ^ 

2 

Y 10 Y j at z ^ z -- 600 (14) 

v,max v,d n 

At redshifts < z , the fluid velocity decays (1+z) because F is 
n 

no longer increasing and Y^ decreases rapidly with z. The decay of the 
fluid velocity and the growth of d (1+z) ^ for z < z^ are due to the cos- 
mological expansion (Ozernoi and Chibisov 1970), The annihilation rate then 

Y nv^(M /M. )d ^ (13) 

V f' a t 

From (A2-19) and (A2-20) , it follows that 

4(v^/c)^ s ^ 4(v^/c) 

so that we find from equations (14) and (15) 


(16) 
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oc where O^Hrsl. (17) 

The constant of proportionality is determined from equations (9) and 
(14) at z ^ z . The resultant values for v_,d, and Y as a function of z 
are shown in figures 2, 3, and 4 respectively, 

F;pom the upper and lower limits on d , we find that at the end of the 
plasma period, the emulsion regions of pure matter and pure antimatter con- 
tain masses of the order of 

~ 10^^ s (M/M^ s ~ 6 X 10^^ (18) 

33 

where 2 X 10 g is equal to one solar mass. 

Thus, even assuming that only 5-10 per cent of the total mass of the 

universe goes into galaxy formation as is consistent with the turbulence picture 

and the y-ray results (see references in appendix 1), it appears likely 

that the emulsion regions will grow large enough by the end of the plasma 

period to contain masses of the size of galaxy clusters - particularly noting 

3 

the uncertainty in d ( M d ) and the fact that our estimate of d is an 
average emulsion size so that somewhat larger sizes can be expected to exist. 
Thus, in this picture, one would expect that all galaxies in a cluster would be 
of the same type, i.e., all matter or all antimatter. 
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IV: THE ONSET OF TURBULENCE AND GALAXY FORMATION 

The size scale of the emulsion regions and the magnitude of the fluid 
velocities estimated in the previous section to be induced by annihila- 
tion pressure are of sufficient size to result in large scale turbulence. 
Turbulence can be sustained if the Reynolds number of the cosmic fluid on 
a scale I 

R = v^>t/v (19) 

2 

is of the order of 10-10 (Heisenberg 1947) where v is the viscosity of 
the fluid. 

During the plasma period, the plasma and blackbody radiation field 
are coupled and behave as a single fluid with the fluid viscosity determined 
by the radiation fiedd. The coefficient of viscosity is then given by 

4 cem 

\ ^ 15 a^p(p 

(see,e.g, ,Silk 1971) . 

During the plasma period, c = p (l+z)/(l+z ) < p and 

eq 

V ^ ~ (c/p) 10^^(l+z) ^ cm^s ^ (21) 

r ISa^n *^ 

Using as given by equation (21) and the values for v^ and d est- 
imated in the previous section to determine the value of the Reynolds num- 
ber for the fluid, we find that turbulence on a scale O.ld ^ ^ < d can 

2 3 

set in during the plasma period for ^ 5 X 10 z ^ 3 X 10 .We thus would 
expect that large scale turbulence existed 'in the cosmic fluid by the end 
of the plasma period. During this period, the turbulence is subsonic with 
the velocity of sound in the fluid being determined by the radiation field 
to be of the order of 10^^ cm/s (see appendix 3). During the plasma period 
then v^ < u^^ and the plasma behaves as an incompressible fluid. 
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The Kolmogoroff law for incompressible isotropic turbulence has been 
well verified in the laboratory and has been nicely summarized by von Weiz- 
sScker (1951) for astrophysral applications. We will briefly summarize here 
von WeizsScker ^s discussion. 

The energy density dissapated per unit density per unit time T through 
the effect of fluid viscosity is given by 

Fa I V X (22) 

where v^ is the turbulent velocity of an eddy of the fluid. 

An eddy of size t has an effective eddy viscosity 


'■‘i - ''i* 


From dimensional arguments, it follows that 


V X vpj} 


so that from equations (22) through (24) , we find 


r - 


(23) 


(24) 


(25) 


Since, in equilibriim, the energy flow through eddies of all sizes 
must be a constant, it follows that 


and 




1/3 




4/3 


(26) 

(27) 


The minimum eddy scale below which turbulence is rapidly dissipated 
by viscosity is that for which the Reynolds number is ~ 1 , so that from 
equations (19) and (25) we find 


( 28 ) 
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The Kolmogoroff spectrum given by equation (26) is established for all 
eddy scales 

^ I . ^ I ^ ^ d (29) 

min 

At the end of the plasma period (z =* z^) when neutralization occurs , 
the radiation field decouples from the matter and two drastic changes occur: 
(1) The viscosity of the fluid drops dramatically to the kinematic viscosity 


10"^^ ^ 3 X 10^"^ (1+z)'^ cm^/s (30) 

(Spitzer 1962, Zel'dovich, et al. 1968). This drop in the viscosity results 
in a very large Reynolds number for the large scale fluid motions and exten- 
sive turbulence down to scales t ^ d. (2) The velocity of sound in the 

min 

fluid drops by about four orders of magnitude (see appendix 3) so that v^ > u^ 
and the turbulence becomes supersonic. The fluid pressures resulting are then 
greater than the kinetic gas pressure and the fluid becomes compressible. Real 
shock waves of large dimensions cannot exist, however, because density fluctu- 
ations on the smaller eddy scales immediatdy divert the wave fronts into stat- 
istically distruibuted directions . (von Weizs^cker 1951). The net result of the 
supersonic turbulence is to cause the large scale density fluctuations needed 
to trigger galaxy formation in an expanding universe (von WeizsScker 1951, 
Gamow 1954, Ozernoi and Chibisov 1970). 

At the end of the plasma period (z ^ 2 ^), as we showed in the previous 
section, the maximum fluid velocity on a scale d can be estimated to be 

1.5 X 10^ cm/s< V (z )< ~ 1.5 X 10^ cm/s (31) 
max n 

During the neutral period, z < z^ , Kolmogoroff 's law for an in- 
compressible fluid (equation (26)) must be modified to allow for compress- 
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ibility and becomes 

where, experimentally 0 < p, < ~ 0, 12 (33) 

and a good estimate for |j, under astrophystal conditions is 

[i^O.07 (33) 

(von WeizsScker 1951). 

The maximum turbulent velocity as a function of redshift will not 
decrease faster than v cc (l+z) (Ozernoi and Chibisov 1970). Thus 


where 1 ^ C js ^ 10 and v 

n 


1 0 ^ 

Vj(-t,z) = Cv^ (L/D) ■ (35) 

n 

g 

1.5 X 10 cm/s, and the comoving variables 


L = >f(l+z) 
D = d(l+z) 


(36) 


According to Ozernoi and Chibisov, gravitationally bound protoclouds 

can form when the virial theorem can be satisfied. The turbulent kinetic 

2 

energy per gram decays qc (1+z) and the graviJational potential energy per 
gram scales as oc (l+z). Ozernoi and Chibisov show that (2T) + ( U ) < 0 
for z ^ z- where 

D 


At that point, the maximum velocity 

~ 2 X 10^ cm/s s V (z. ) s ~ 2 X 10^ cm/s (38) 

max b ^ 

«3 

and the gas density n(z^) 1 cm so that 1 galactic mass is contained 

within a scale 

I =“ 3 X 10^^ cm 
g 


( 39 ) 
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for which the maximum rotational turbulent velocity from equation (35) is 


2 X 10^ cm/s s V (^ ) 5 ~ 6 X 10^ cm/s 
max' g' 


(40) 


After , the protogalactic clouds maintain their integrity as in- 


dicated by the fact that for our galaxy 
and the galactic radius 


n - n(z^) 


r ~ -C, 
g S 


(41) 


(42), 


Further galactic evolution may follow along the lines suggested by von Weiz- 
s^cker (1951) but will not be examined here. 

We also find observationally that galaxies have a maximum observed 
rotational velocity 

V , ^ 3 X 10^ cm/s (43) 

rot ,max ' 

(Brosche 1967) in agreement with equation (40). 
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V: THE y-RAY BACKGROUND SPECTRUM: 


The most direct empirical verification of the extension of the Omnes 
cosmology suggested here lies in a study of the cosmic y-ray background 
spectrum above 1 MeV ; those y-rays being presumably the direct products of 
the decay of rr^-mesons produced in proton-antiproton annihilations (Stacker, 
et al. 1971). Two characteristics of the spectrum have been determined by 
Stacker, et al.: (1) a flattening of the y-ray spectrum in the neighborhood 
of 1 MeV due to absorption effects at z ^ - 100 and (2) a power-law form 
for the y-ray spectrum at higher energies ( ~ 5-50 MeV ). 

If we denote the unredshifted y-ray spectrum fron proton-antiproton 


annihilation by 


then the y-ray spectrum predicted for an Einstein- 


de Sitter universe (see appendix 1) is given by 


(1+2)“^’^ t(z) f^[E^(l+z)] (44) 

(Stecker 1971), where the integral is taken over redshifts where absorption 
is unimportant (z^ ^--^100) so that equation (44) is only valid for ^ 5 ^ E < 

y 

~ 50 MeV (Stecker, et al. 1971), 

Ifj we then represent Y(z) in the form 

f(z) oc ( 1 +z)^^ (45) 

as suggested by equation (17), then 


where 


I^(E^) «Jdz (l+z)“ f^^^(l+z)] 
(X = 1.5 + H 


( 46 ) 


Then, because of the bounded form of f^(E 



the cosmological y-ray 


spectrum between ^ and ^--50 MeV has the power -law form 
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I. (E ) dE oc E dE 

A Y Y Y Y 

- (2 5+h) 

oc E jjg 

Y Y 

where 0 ^ k ^ 1 

(Stecker 1971), 

For the microscopic model discussed by Stecker, et al. (1971), the 
comparable value of h was 0.36 ; the macroscopic model gives values for h 
between 0 and 1. In the notation of this paper, the effective annihilation 
rate used by Stecker, et al. which was normalized to the y-vsy observations 
IS (in cm s ) 

Y = 0.8 X 10"^'^ (l+z)^’^^ z ^ ~ Zu (48) 

Y u 

which, of course, is only valid for the lower redshifts z ^ z^. where the 
universe is transparent to y“i^^di^tion. An extrapolation of this function 
to higher redshifts leads to rotational velocities ~ 10^ cm/s for our turb- 
ulence model which compares well with the observed value as given by equa- 
tion (43). It is also within order of magnitude agreement with equation (9). 
Gravitational binding and subsequent strengthening of magnetic fields for 
z ^ z^ may be expected to reduce the annihilation rate somewhat below the 
extrapolated value flfom higher redshifts. 

We conclude that the y-ray evidence as discussed by Stecker, et al. 
lends support to the picture presented here for galaxy formation by 
annihilation -generated turbulence. Other effects of the model presented 
here may be (1) some distortion of the 2.7 K blackbody spectrum on the 
high-frequency side, and (2) heating of the intergalactic medium by the 
decay of turbulence at low redshifts. 
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Appendix 1: Cosmological Relations Used in the Calculation, 

We adapt here for the geometry of our cosmological model, the isotropic 
and homogeneous solution to the Einstien field equations givenfor the Robertson- 
Walker metric 

7 2 2 2 2 

ds = c dt - R""(t)du , (Al-1) 

where R(t) is the scale with which the universe expands as a function of time t 
after 'the big bang. If a photon is emitted at a time t and received at a time t^, 
it can easily be shown that the relation between the amount the photon has been 
Doppler-shifted due to the expansion of the universe, AX/X ^ z, and the time of 
emission is given by 

R^/RCt) = 1 + z. (Al-2) 

It then follows from (Al-2) that particle densities at time t scale with redshift 
z as follows: 

Photon energies scale with redshift as the frequencies according to the Planck 

relation E = hv so that 
Y 

E = E (1 + z). (Al-4) 

Y Y^O 

From (Al-3) and (Al-4), it can be shown that the temperature of the blackbody 

radiation scales with redshift as 

T - Tq(1+z) (Al-5) 

where is observed to be 2.7 K. It follows from (Al-5) that the photon energy 

4 4 

density e cc X « (1+z) . 

Y 

The homogeneous -isotropic solution to the Bins tan field equations contains 
only two arbitrary constants which must be determined by observation (if, as here, 
we set the"cosmological constant” x = 0) ; they are the Hubble constant, , and 
a parameter Q which is proportional to the mean matter density in the universe. 
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We will choose here for the Hubble constant, the most recent value of 
~ 1.7 X lO” sec (Abell 1971), The parameter Q is a measure of the deceleration 
of the universe due to gravitational effects and is defined as 

Q = (A 1-6) 

where n is the mean density of matter needed to close the universe gravitationally, 
c 

In this discussion, we will take Q = 1 in order to preserve consistency with the 
antimat ter -annihilation explanation of the observations (Stacker, et al, 

1971) and with the turbulence hypothesis for galaxy formation (Oort 1969). This 
corresponds to the Einstein-de Sitter model (which, in any case, is valid for Ozl>l) 
which gives the time-redshif t relation 

t = tQ(l+z)“^/2 (Al-7) 

where t^ is the present age of the universe (<^ 2 X lO^^yr). 

A more complete discussion may be found in references such as McVittie (1965). 
Simplified derivations for application to y-ray astronomy may be found in Stacker 
(1971). 


Appendix 2: Annihilation of Matter and Antimatter in Boundary Regions. 


A, Annihilation rate when the fluid velocity v^ is less that the thermal 
velocity v^ . 

In this case, mixing of matter and antimatter by diffusion determines the 
mean annihilation rate. The equations for matter and antimatter density and 
N- as a function of distance x paerpendicular to the boundary surface are 


dN 

^ = 

- < 


> N N- + D 
' P P 

a^N 

dN 

dN- = 
at^ 

- < 

^A^t 

> N N- + D 
^ P P 

a^N- 
^ - 

aN- 

'^d^P 


(A2-1) 


where a is the cross-section for matter -antimat ter annihilation, D is the 
A ’ 


effective diffusion coefficent and v^ if the effective diffusion velocity 


for matter and antimatter erlEring the annihilation region (Bessard, Dennefeld, 
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and Puget 1969). The stationary solutions of (A2-1) are given by 


2x/h 


N s 
P 0 


N- = N 
P 0 


4 cosh (x/h) 

-2x /h 
e 

4 cosh^(x/h) 


(A2-2) 


with the boundary conditions N^(oo) m N- (-<») ^ and = N-(<») - 0 (Omnes 

1971^) .The characteristic dimension h, and the diffusion velocity v^ are given by 

h 2(6D/N^<o^v^»i 


and 


(A2-3) 


v^ = 2D/h . 

G 

The total annihilation rate per unit surface area is then 


2 

dx <a.v >N N- ^ ^0 ^ . (A2-4) 

„ A p p 

12 


The diffusion coefficient is given by 


so that 


and 


D = i 


h = 2/2 X, 


Vd ^ v^/V2 


(A2-5) 


(A2-6) 


where is the annihilation mean-free -path , equal to (Na^) 
The resultant expression for then simplifies to 


Yg = 0.37 


(A2-7) 


B. Annihilation rate when the fluid velocity v^ is greater than the 
thermal velocity v^. 

, When v^ > v^ , macroscopic masses of matter and antimatter will collide as 
a result of the random fluid motions and the complicated geometry of the emul- 
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sion and considerably disrupt the stationary state just described. This situation 

4 

occurs for z < 4 X 10 (see figure 2). Let us consider such a collision beginning 
at time t 0. For ^ ^ matter and antimatter mix on the boundary 

with a rate Nv^ per unit area. For times t , a strong flux of y rays and rel- 

ativistic electrons is emitted at the boundary region. We consider here the energy 
carried by the electrons which have an average energy of -- 100 MeV and neglect the 
effect of the y-rays whose mean f ree -^path is of the order of the emulsion size, d. 
The electrons lose their energy through Compton collisions with blackbody photons 


and produce X-rays with an average energy 


30TT2)2 (2*7kT-) - 3S(l+z) eV (A2-8) 


The lifetime of the electrons is then 


(A2-9) 

where n is the photon density of the blackbody radiation and is the 
Thompson cross section. 

The X-rays in turn lose their energy through Compton collisions with 
electrons of the cosmic plasma and their mean-free -path is 

= (n = 5 X 10^^ (l+z)"^cm . (A2-10) 

A P 

In each collision, the X-rays lose.:an energy given roughly by 

AE ssE^/m sa 2 X 10"^(l+z)^ eV . (A2-11) 

The time for transport of momentum from the X-rays to the plasma is 

then 

T X^/c ^ 1.7 X 10^^(l+z)“^ s (A2-12) 

mom X 

and the time for transport of energy is given by 

. Tg (X^/c)(E^/AE^) « 2.7 X 10^^(l+z)'^ s. (A2-13) 
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The average momentum transmitted per second during the burst of annihila- 
tion is 


^ -27111 c/n ) a 0.27m V cX' 

At ^ P P P^ ^ (A2 14) 


where the factor 0.27 is the fraction of annihilation energy transmitted to the 
electrons and positrons which are among the final products of the annihilation 
of the protons and antiprotons. The annihilation rate per unit area is nv^ for 
times t > T* because t > t* for most of the period we consider. The change 
in the velocity of a nonrelativistic electron due to electrostatic interactions 
with protons is given by 






(A2-15) 


(Spitzer 1962) where v^ is the electron velocity, v^ is the thermal velocity of 
the protons , ^n ^ and G are functions of v^, v^ , and n. From (A2-14) and (A2- 
15) we see that the electrons and protons must have a relative velocity v^ given 
by the relation 


0.17 V/f WpC 


- 4 


(^A) &( ^ ) 


(A2-16) 


5 4 

Eqaation (A2-16) gives v^ ~ 10 cm/s for z ~ 10 . Such a relative velocity cor-- 
responds to a hugh electric current and will rapidly lead to fields which will 
stop the relative motions of the electrons and protons and transmit the momentum 
of the electrons to the whole fluid. With such fields , ^6 o? Ap/At from equation 
(A2-14) and fluid motions with velocity v^ toward the boundary will stop in a time 


such that 

which implies that 


Set - m v - 
P f 


t 


X^/0.27c 


4t 

mom 


(A2-17) 

(A2-18) 
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Thus, if we consider a mass of fluid moving with velocity v^ , it will be 

stopped by annihilation pressure at the boundary region in a time comparable to 

T within a distance \ of the boundary with the annihilated mass M 'Such that 
mom X a 

the momentum produced is equal to the initial momentum of the gas. Denoting the 
total mass of the fluid within a distance \ by M , 


0.27 M c = M V _ . 
a X f 


(A2-19) 


For fluid motions involving total masses M > M , one can still consider the 

^ X 

mass of the fluid within X^ show that equation (A2-19) still holds. 

For more complex fluid geometries than we have considered here, in general 
d^/St ^ 0 and magnetic fields will be created such that they will be amplified 
by the fluid motions to reach a large enough value to cushion them at the bound- 
ary (Schatzman 1970, Puget, to be published). When the magnetic field plays a 
dominant role, the condition of equipartition between the energy of the field 
produced and the energy going into heating the plasma leads to the relation 


%(0.27M^c^) = (A2-20) 

which gives a mueh smaller value for M than equation (A2-19). Thus, most of the 

a 

annihilations will occur in interactions when the magnetic field does not play a 
dominant role. 

We now consider the transfer of energy from the annihilation into macro- 
scopic fluid motions. Electrons in the boundary region are heated at a rate 


(dE^dt) = (0.27iripC^v^/0.27c) (A2-21) 

as obtained from equation (A2-19) where is given by equation (A2-13). 

Assuming the emulsion size d scales like X (see main paper and Omnes references) 
so that V- scales like dX /dt) . 

The y-ray mean-free-path 

\ s (n(j ) 1.6X10 (1+z) cm 

It then follows from equations (A2-13) , (A2-21) and (A2-22) that 


(A2-22) 
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(dE^/dt) =^9.1 X 10”^^ (1+z)^^^ eV/s (A2-23) 

X 

The energy loss rate for cooling of electrons by interactions with the 

blackbody radiation is given by 

(dE^/dt)j^l^ =“ ,8/3 a^c (E^/m^c^)e (A2-24) 

-90 4 

- 2.7 X 10 (1+z) eV/s 

Weymann (1965). 

It follows from equations (A2-23) and (A2-24), it follows that antimatter 
annihilation will heat up the electrons in the cosmic plasma at the boundary 
regions to an equilibriiim energy given by 

E ^ 3.4 X 10^ (1+z)"^^^ eV (A2-25) 

eq 

At this electron energy, energy is being fed through the plasma in a 
characteristic time obtained from equation (A2-19) 

T = m V c (dE /dt ) ^ 

c p f G hb 

(A2-26) 

We also note the other condition on the heating of the plasma electrons, 
viz., that they can only be heated to a maximum energy by the X-rays as given 
by equation (A2-8). The allowed energy range for the plasma electrons as given 
by the restrictions (A2-8) and (A2-25) is shown in figure 5 as a function of 
redshift. It follows from these restrictions, as shown in figure 3 that cooling 
by the blackbody radiation is unimportant for redshifts 600 because in this 

redshift range, electrons cannot be heated to E as given by equation (A2-25). 

The hot gas of electrons transmits its energy to the protons within a 
time given by (Spitzer 1962) 

Y S (A2-27) 

P 32tt M A Cr( 


1. 1 X 


S 
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where w is the velocity of the hot electrons and -tnA.and G(w /w ) are the same 

6 © p 

functions as in (A2-6), It follows from equation (A2-27) that ^ 2 X 10^^ s 


for (1+z) ^ ~600. 


The hot gas can expand almost freely within a dynamical time scale t 
given by 


dyn 


dyn 




(A2-28) 


where is the thermal velocity of the hot protons in equilibrium with the hot 

electron gas. If t < t and t < t, , condition which is fulfil- 

^ pc p dyn 

led for all z ^ 600, then E ^ E and 

P e 


Vp = (2E^/mp)^ =-2.6 X 10^°(l+z) (A2-29) 

The dynamical expansion of the hot gas produces macroscopic fluid motion 
in the directions away from the boundary. Thus, a fraction F = 

energy released in the form of relativistic electrons in the annihilation is 

transformed into macroscopic fluid motion rather than being cooled by the 

blackbody radiation field. From equations (A2 -24) , (A2-26) , (A2-28) and (A2-29), 

we then obtain an expression for the dynamic efficiency F 


F 7 X 10^ (1+z)'^'^'^ for (1+z) & ~ 600 (A2-30) 

and F ==* 1 for (1+z) ^ ~ 600 (A2-31) 

when cooling by blackbody radiation is no longer important. 
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A p p e nd lx 3 : Ne u t'r a 1' i z a t ion Per rod of * the C ds m i c P las ina ; ' ^ ‘ ■ 

The neutralization of the cosmic plasma is slower in the baryon -symmetric 
cosmology than in the normal "big-baoig" cosmology -because the annihilation prod- 
ucts keep the plasma ionized at a higher level than would be expected from pure- 
ly thermal considerations, A detailed solution of the neutralization problem 
for the normal "big -bang" cosmology has been given by Peebles (1968). The equa- 
tion he given for the electron density is 

dOe/dt = C(n^,T) C exp -(Bj^-B^)/kT } (A3-l^ 

where aBO ‘ ■ is the b ind ing energy for ; a;>hydr ogen a tom , with i t s ! e 1 ec t r oh - in an orb it / 
of quantum number n and n^ is the density -of hydrogen atoms in the ground state. 
The and 3^ coefficients, are functions :of ‘ temperature T, - " 

In order to investigate the effect of the annihilation products on the 
neutralization of the cosmic plasma, we have considered the behavior of these 
products. The neutrinos which carry roughly half of the annihilation energy 
have no ionizing effect. The lose their energy by pair production and 

Compton interactions with an effective cross section ^ 2 X 10 cm (See, e.g.. 
Stacker 1971) with the secondary electrons having an energy of the. same order 
of magnitude as that of the parent - y-rays are also produced by the 

annihilation of electrons and positrons after they have been slowed down (see 
also Stecker 1971). The electrons and positrons resulting from matter antimatter 
annihilation, of course, produce X-rays as we have discussed in appendix 2. 

The X-rays can themselves undergo two interaction processes, viz., 
photoionization and Compton interactions. The photonionization cross section 
depends strongly on and is given by 


(A3 -2) 
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where cohere is the fine structure constant . (see, e.g. , Heitler 1954). 

The X-rays are very efficient ionizing agents and tend to slow down the 
neutralization of the plasma. At the beginning of their life, the most energetic 
X-rays are degraded by Compton interactions with free electrons. This process, 
however, becomes unimportant when the energy of the X-rays falls below a few keV, 
at which point photoionization becomes important if only 1% of the hydrogen is 
neutral. 

The number of ionizations per annihilation is thus given by 

=>■ 0.5 (2mpC^/32eV) (1-x) (A3-3) 

(where x is the fraction of hydrogen in an ionized state) because one ionization 
corresponds to the ultimate loss of 32 eV by a fast electron. 

The equation for ionization equilibrium then becomes 

dn^/dt = a^nl CCn^^.T) - k. a 0 (A3-4) 

where ¥ is the annihilation rate per unit volume. 

Equation (A3 -4) reduces to 

C(T,x)x^ = 6.75 X 10"^ (1-x) (1+z) (A3-5) 

if we adopt the value given by equation ( 9 ) of the main text as a lower limit 
3 

on at z ^ 10 (see main text for discussion). 

The solution to equation (A3 -5) is shown in table A3-1. It follows from 

this table that “total" neutralization takes place later in the baryon- symmetric 

cosmology than in the normal all*matter cosmology. Several remarks are in order 

3 it 

concerning this result. For T < 10 K, x ocx® . However, this value of x is an 

average; stronger ionization will occur in the vicinity of the boundary regions. 

Another important point is that the annihilation rate decreases rapidly with z. 

As we show in the main text, part of the annihilation energy may be stored in 

3 

large scale turbulence for 500 < z < 10 . Part of this kinetic energy is therm- 
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alized later when the recombination rate ( proportional to n ) is much smaller. 
Thus, we should consider the results of table A3-1 to represent lower limits on 
the value jof x. 

The sound velocity for a plasma coupled to the blackbody radiation field 




(A3 -6) 


where n is the blackbody photon density and n is the proton density. 
Y P f 

The sound velocity for a neutral gas at temperature T is 


u a (5kT/3m ^ (A3 -7) 

n ' p' t 

The ratio 

u /u , ~ (n /n )8 ^ (A3-8) 

n pi P Y 

Thus, the sound velocity of the cosmic gas drops by about four orders 
of magnitude by the end of the neutralization period and the fluid velocities 
considered in part III of the main text will then become supersonic (see 
figure 2 ), The transition takes place when the mean- free -path of the photons 
becomes larger than the characteristic fluid dimension d.(see, e.g., Ozernoi 
and Chernin 1968 for further discussion). The mean free-path of the thermal 
photons is 

^th “ (A-3-9) 

where X is given by equation (A2-10). The X-rays, as shown in appendix 2, have 
a lifetime much longer than X /c.If this lifetime is^ the hydrodynamic time t, = 
d/v^ , they will contribute to maintaining a high sound velocity (The effect of 
high-energy photons on the sound velocity has been studied by Montmerle 1971). 
This condition for the high sound velocity is fulfilled for z ^ ~ 500, Figure 
3 shows that the mean -free -path for the thermal photons becomes larger than 
d for z ^ ~ 600, We will thus consider neutralization (corresponding to the 

transition from u^^ to u^ to occur for z^ 5±1 X 10^ (keeping in mind that we 
may have somewhat underestimated the value of x(z). 



TABLE 3A-1 . NEUTRALIZATION OF THE COSMIC PLASMA 


z 

Peebles (1968) 

This Work 


X 

C(T,x) 

X 

C(T,x) 

1660 

0.92 

i.8 X lO”^ 

0.96 

3.6 X 10“^ 

1480 

0.40 

5.9 X 10"^ 

0.84 

2.2 X 10"^ 

1300 

0.072 

2.7 X 10"^ 

0.64 

7 X 10"^ 

1100 

9.8 X 10“^ 

2 X 10"^ 

0.38 

3.2 X 10"^ 

925 

9.2 X 10“^ 

0.-25 

0.14 

0.3 

740 

1.23 X 10"^ 

0.96 

! 0.068 

1 


550 


5.3 X 10 


0.06 


370 


0.049 
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FIGURE CAPTIONS 


Figure 1* The growth of the matter and antimatter regions as a function 
of cosmic temperature ( T = 2.7K(l+z)). 

Figure 2. Various velocities discussed in the text given as a function 



plasma , spped of sound in the neutralized gas, v^ as given 

in equation (12), v^ as given in equation (10) which gives the lower 
limit on the fluid velocity at z^ ; v^ is the thermal velocity and 
v^aj^(l) given by equation (1). 

Figure 3, Range of emulsion size d as a function of redshift as discussed 
in the text. 

Figure 4. Various annihilation rates discussed in the text as a function 
of redshift. 

Figure 5. Allowed range of electron energy as determined by equations (A2-8) 
and (A2-25). 
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